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ABSTRACT 

Infrared emission from the dust shell around IRC+10216 is analysed in detail, em- 
ploying a self-consistent model for radiatively driven winds around late-type stars 
that couples the equations of motion and radiative transfer in the dust. The resulting 
model provides agreement with the wealth of available data, including the spectral 
energy distribution in the range 0.5-1000 |^m, and visibility and array observations. 
Previous conclusions about two dust shells, derived from modelling the data with a 
few single-temperature components of different radii, are not supported by our re- 
sults. The extended, continuous temperature and density distributions derived from 
our model obviate the need for such discrete shells. The IR properties vary with the 
stellar phase, reflecting changes in both the dust condensation radius ri and the overall 
optical depth r - as the luminosity increases from minimum to maximum, ri increases 
while T decreases. We find that the angular size of the dust condensation zone varies 
from 0.3 arcsec at minimum light to 0.5 arcsec at maximum. The shortage of flux at 
short wavelengths encountered in previous studies is resolved by employing a grain size 
distribution that includes grains larger than 0.1 |^m, required also for the visibility 
fits. This distribution is in agreement with the one recently proposed by Jura in a 
study that probed the outer regions of the envelope. Since our constraints on the size 
distribution mostly reflect the envelope's inner regions, the agreement of these inde- 
pendent studies is evidence against signiflcant changes in grain sizes through effects 
like sputtering or grain growth after the initial formation at the dust condensation 
zone. 

Key words: stars: late-type - stars: mass-loss - stars: individual (IRC+10216) - 
infrared: stars ~ circumstellar matter - dust 



1 INTRODUCTION 

IRC+10216 (CW Leo, IRAS 09452+1330) is by far the 
brightest and best-studied mass-losing carbon star (Jura 
& Kleinmann 1989). Starting with Mitchell & Robinson 
(1980), several authors have performed radiative transfer 
calculations for the IR dust emission from this source. How- 
ever, with the exception of Winters, Dominik & Sedlmayr 
(1994), all the previous studies were based on a prescribed 
radial density distribution that is not fully consistent 
with those of radiatively driven winds. Furthermore, al- 
though reasonable fits to the observed spectral energy dis- 
tribution (SED) were generated over a wide range of wave- 
lengths, none of the models produced enough flux shortward 
of 1-2 \im (e.g. Le Bertre 1987; Keady, Hall & Ridgway 
1988; Griffin 1990; Lorenz-Martins & Lefevre 1993). In ad- 
dition, none of the models provides simultaneous agreement 
with spatially resolved observations at 2.2 ^.m (e.g. Martin 
& Rogers 1987). 

The purpose of this work is to perform a self-consistent 



study that employs a dust density distribution determined 
from the solution of the coupled system of radiative trans- 
fer and hydrodynamics equations for the wind. The equa- 
tions are described elsewhere (Netzer & Elitzur 1993; Ivezic 
& Elitzur 1995, hereafter IE95). As shown in IE95, the so- 
lution of this system is essentially determined by a single 
quantity - the flux-averaged optical depth rp. Once rp is 
determined, scaling relations listed in IE95 and in Ivezic & 
Elitzur (1996; hereafter IE96) can be used to constrain all 
other relevant quantities. In principle, the optical depth can 
be estimated from either the spectral shape f\ = F\/F, 
where F = J -FAdA is the bolometric flux, or spatially re- 
solved observations. However, since the latter depend also 
on the angular scale of the system, because of observational 
uncertainties the determination of optical depth from the 
spectral shape is much more reliable. In IE96 we describe 
a two-step modelling procedure, which we follow in this 
work. In the first step the dust characteristics and overall 
optical depth are constrained from the best fit to the spec- 
tral shape. Then, with the model prediction for the surface 
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Figure 1. Spectral energy distribution for IRC +10216; lines 
represent model results, symbols the observations. Data are from 
Le Bertre (1987) (■), (1988) (o); Rengarajan et al. (1985) (•); 
and IRAS Point Source Catalogue (*). All observations are at 
maximum light except for those denoted by open circles, which 
were at minimum light. The thick solid line is the model result 
for maximum light, the thin solid line the result for minimum; 
details are described in the text. The dashed line is the model 
result for maximum light and single-size (0.05 ytsn) grains. The 
inset shows an expanded view of the IRAS LRS spectral region 
— the dots are the data, taken close to maximum light, the solid 
line the model. 

brightness distribution based on these parameters, the spa- 
tially resolved observations are used in the second step to 
determine the angular size of the dust condensation zone. 

The spectral energy distribution is discussed in Section 
2, high-resolution observations in Section 3 and outflow dy- 
namics in Section 4. The results are summarized and dis- 
cussed in Section 5. 



2 SPECTRAL ENERGY DISTRIBUTION 

IRC-l-10216 is a long-period variable with a period of 638 
d and a recent minimum at JD = 244 7863 (Dyck et al. 
1991). We limit our analysis to the periodical changes and 
do not consider the long-term modulations with time-scales 
of a few decades noted by Dyck et al. The surface bright- 
ness distribution at short wavelengths (sC 2-3 \im) is slightly 
asymmetric (e.g. Ridgway & Keady 1988; Kastner & Wein- 
traub 1994). However, these spatial asymmetries decrease 
as the wavelength increases (e.g. de Batz 1988)[| and our 
model assumes spherical symmetry. Slight disagreement be- 
tween the model and observations can be expected at short 
wavelengths, a point further discussed in Section 5. 

Our best-fitting model to the spectral shape is shown 
in Fig. 1 together with the observations. The thick solid 
line corresponds to maximum hght, the thin solid line to 
minimum (where there are only five observational points). 
The inset compares the model results (solid line) with the 
IRAS LRS data. The model is primarily determined by the 
overall optical depth and the dust composition.^ From pre- 
vious work (e.g. Blanco et al. 1994), the dust grains around 
IRC-l-10216 are primarily composed of amorphous carbon 
with a minor inclusion of SiC to account for the 11.3-ii.m 
feature. With optical properties for amorphous carbon taken 
from Hanner (1988) and for SiC from Pegourie (1988), we 
find that the best fit to the 11.3-[j,m feature is obtained with 

* Recently, Sloan & Egan (1995) observed IRC-l-10216 with a 
0.9 X 2.0 arcsec^ resolution at 10 \Ma and obtained an indication of 
a blue emission patch with a size of ~ 1 arcsec, located ~ 1 arcsec 
north of the star. Note that these spatial details are comparable to 
the slit size; indeed, Sloan & Egan point out that this component's 
location could be an artefact of the reconstruction algorithm and 
that it could actually originate from the inner region centered on 
the star. Such emission indeed is expected from this region due 
to the hot dust. These observations do not seem to refute the 
findings of de Batz. 

t The required grain properties are the spectral shapes of the 
absorption and scattering efficiencies. Absolute values of these 
quantities are not needed. 



a mixture of 95 per cent amorphous carbon and 5 per cent 
SiC (by mass), although varying the percentage of SiC in 
the range 3-8 per cent still produces satisfactory agreement. 
Griffin (1990) presents results for various SiC abundances. 
He obtains the best agreement for 17 per cent SiC, but the 
quality of his fit at 8 per cent is comparable. The slight 
difference between the conclusions of the two studies is in- 
significant and might be explained by the density law 
employed by Griffin instead of the hydrodynamic calculation 
done here. A broad emission feature between 24 and 30 \Lm 
(Forrest, Houck & McCarthy 1981) provides evidence for an 
additional component, probably MgS compound (Goebel & 
Moseley 1980). With the aid of spectroscopic data for MgS 
from Nuth et al. (1985) we estimate the abundance of this 
component to be less than 10 per cent, if this chemical iden- 
tification is correct. 

In addition to the chemical composition, the distribu- 
tion of grain radii a also affects the optical properties. How- 
ever, the wavelength dependence of absorption and scatter- 
ing efficiencies is independent of a once A ^ 2na. Therefore, 
at the wavelengths of interest, A 0.5 \im, the grain size is 
irrelevant as long as a <, 0.1 [im, and models of IRC-l-10216 
usually assumed that all grains have the same size a — 0.05 
\im (for an overview see Lorenz-Martins & Lefevre 1993). 
The dashed line in Fig. 1 shows our model result for this 
single size, displaying the problems encountered by all other 
workers - the models do not produce enough flux at A ^ 1 
\im. 

What is the meaning of this discrepancy? Dust emis- 
sion is insignificant at A ^ 3 \im because it decreases expo- 
nentially for wavelengths shorter than 3 [im x (1000 K/Ti), 
the peak wavelength of the Planck distribution for the dust 
condensation temperature Ti , the highest possible dust tem- 
perature. Therefore, the detected radiation involves only at- 
tenuated stellar emission and scattered light. Our detailed 
models show that, under these conditions, the spectral shape 
is proportional to exp(— Tabs), where Tabs is the overall op- 
tical depth for absorption. Therefore, the shortage of ob- 
served fiux at A 1 [i-m implies that the model estimates 
for Tabs (A) at these wavelengths are too large. Since in gen- 
eral Tabs (A) increases as the wavelength decreases, this rise 
must be suppressed around 1 \Lm. Indeed, Rowan- Robinson 
& Harris (1983) noted that by postulating a departure from 
the usual A"'^ dependence of rabs(A) to a fiatter distribution 
at A ^ 1 ^m they could produce a better agreement with 
observations. Although they did not attempt to justify this 
behaviour, it can be modelled by assuming a range of grain 
sizes a. While Qabs oc A^^ for A <; 27ta, it is approximately 
constant for A < 27ta. Consequently, the fiux shortage at 
A ^ 1 |i.m can be alleviated by adding larger grains with 
sizes of a ~ A/27t ~ 0.2 (i.m, suppressing the rise of rabs(A) 
when A decreases below ~ 1 ^tm. 

For the detailed models we employed two types of size 
distributions n{a). Most often used is 

(1) 

proposed by Mathis, Rumpl & Nordsieck (1977). This MRN 
distribution includes a sharp cutoff, flmax, to the grain radii, 
required by the finite amount of mass in the dust. Recently 
Jura (1994) proposed a modification of the form 

nj(a) oc a-^^-^e""/"", (2) 
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replacing the sharp cutoff with an exponential one. Both dis- 
tributions can produce satisfactory fits: the MRN distribu- 
tion requires Omax ~ 0.2-0.3 \i.ra, the Jura distribution ao ~ 
0.15-0.2 ^.m. It appears that the Jura distribution produces 
a slightly better fit to the spectral shape observed at maxi- 
mum light, but any stronger conclusion is hampered by the 
observational uncertainties. In contrast, recent models by 
Bagnulo, Doyle & Griffin (1995) produced satisfactory fits 
with single-size grains of 0.02 fi.m as well as the MRN distri- 
bution with flmax = 0.05 li-m, but not with the Jura distribu- 
tion. However, these models used an density law instead 
of employing a self-consistent hydrodynamic calculation, as 
done here. The two density distributions are substantially 
different at the inner regions - the self-consistent distribu- 
tion has a much faster initial fall-off and has already dropped 
by a factor of 5 below the distribution at r ~ 1.5ri . Since 
this is the region where the short wavelengths are produced, 
this could account for the different results. 

We have thus determined the two major ingredients 
that affect the spectral shape, the grain optical properties 
and overall optical depth. In addition, the stellar tempera- 
ture T* and dust condensation temperature Ti have a dis- 
cernible effect on the spectral shape, but only at short wave- 
lengths. Our best fit gives T* = 2200 ± 150 K, an estimate 
in agreement with a spectral type of C9 (Cohen 1979) and 
the majority of other models. In general, the effect of T* is 
limited to A 4 ^m and its significance is diminished as the 
envelope's optical depth increases. Our best-fitting estimate 
for Ti is 750 ± 50 K. The effect of T\ is more significant 
because this parameter controls the peak wavelength of the 
spectral shape in envelopes that are optically thin around 
that peak. Our estimate for Ti, determined from the ob- 
served spectral shape by the location of the peak and the 
sharp decline toward short wavelengths, is somewhat lower 
than the ~ 1000 K obtained in most other models. Indeed, 
in these models the peak of the spectral shape is shifted 
slightly to the left, resulting in excessive fiux in the 2-7 ^.m 
wavelength range (e.g., Le Bertre 1987; Lorenz-Martins & 
Lefevre 1993). 

Finally, the radius of the envelope's outer edge, rout, 
must be specified for a numerical solution. Because of scal- 
ing, only the relative thickness j/out = ^out /f i is needed. This 
parameter affects only the long-wavelength part of the solu- 
tion, which is afflicted by a number of uncertainties. First, 
the behaviour of the absorption efficiency is quite uncer- 
tain at these wavelengths. Usually modelled by a power law 
Qabs oc A"'', the value of (5 is poorly known, typically taken 
as ~ 1-1.5. Next, the long-wavelength tail of the SED could 
contain a significant contribution from free-free emission 
(Griffin 1990). Fortunately, apart from the long wavelength 
part of the SED, the model results are not very sensitive to 
these uncertainties. We find from fits to the spectral shape 
in the wavelength range 100-1000 \ira that (3 varies from 
1.2 to 1.6 for j/out between 600 and 10 000. An independent 
estimate for j/out can be obtained from the extent of molec- 
ular emission. CO observations by Huggins, Olofsson & Jo- 
hansson (1988) indicate that j/out S5 700, and consequently 
/3 ^ 1.2. Indeed, Jura (1983) suggested that (3 = 1.3, a pro- 
posal supported by Le Bertre (1987). The results presented 
in Fig. 1 are for i/out ~ 700 and /3 — 1.3. For this /?, i/out 
can be increased all the way to 10 000 without a significant 
degradation of the fits. 



Table 1. Overall optical depths for the best-fitting models plot- 
ted in Fig. 1. The last entry lists the flux-averaged optical depth 
Tp. The other properties of the models are; dust composition, 95 
per cent amorphous carbon and 5 per cent SiC (by mass); grain 
size distribution given in equation (2), with ao = 0.2 \im\ dust 
condensation temperature T\ = 750 K; stellar temperature T« = 
2200 K. 



A (|Jm) 




^max (b) 


.^min(c) 


0.55 


0.52 


20 


24 


1.0 


0.44 


14 


17 


2.2 


0.40 


4.7 


5.7 


5.0 


0.13 


1.0 


1.2 


10 


0.026 


0.32 


0.40 


100 


1.0 X 10--* 


0.01 


0.012 






2.0 


2.4 



(a) Albedo 

(b) Total optical depth at maximum light 

(c) Total optical depth at minimum light 

The parameters of our best-fitting model are summa- 
rized in Table 1. Note again that the fit to the spectral 
shape /a is obtained without specifying the absolute size 
of the envelope, mass-loss rate, luminosity or distance to 
the star. For given dust grains, the resulting fx is primar- 
ily determined by the overall optical depth. Furthermore, 
the fiux scale never entered the fitting procedure. Actual 
fluxes are obtained from f\ through simple multiplication 
by the bolometric flux F. Comparison of fluxes from our 
best-fitting model with observations gives F — 2.1 x 10~* 
W m~^ at maximum light, in agreement with Sopka et al. 
(1985). In IE96 we show that the bolometric flux and an- 
gular diameter of the dust condensation zone, 9i — 2ri/D 
where D is the distance to the source, are related via 

^^=°-^^n iO-3Wm-0 [-tt) 

where q is a dimensionless coefficient of order unity char- 
acteristic of the model. This coefficient, determined the- 
oretically from the overall solution, depends primarily on 
the grain optical properties and only slightly on T*, Ti 
and overall optical depth. From our best-fitting model for 
IRC-H0216 we find that a = 1.3 for this source, and there- 
fore at maximum light Q'^'^^ = 0.56 arcsec. With an expected 
bolometric amplitude of 1 mag, 0™'" = 0.35 arcsec at mini- 
mum light. These estimates for the angular scale must agree 
with high-resolution observations. 



3 SPATIALLY RESOLVED OBSERVATIONS 

In interpreting the spatially resolved observations of 
IRC-H0216, the source variability must be taken into ac- 
count. As the luminosity varies during the stellar cycle, the 
envelope temperature varies too. Therefore, as noted already 
by Danchi et al. (1990, 1994), the dust condensation radius 
varies during the stellar cycle and ri scales as LP'^ (IE96). 
Because of the movement of the shell's inner boundary, the 
overall optical depth is expected to vary too, so that maxi- 
mum light has minimum r. With optical depths determined 
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Figure 2. Visibility functions for IRC+10216. Lines represent 
model results, symbols the observations. Solid lines and full sym- 
bols (including + and *) correspond to phases close to maximum 
light, open symbols and dashed lines to phases close to minimum. 
Data are from Sutton, Betz & Storey (1979) (*), Selby, Wade & 
Sanchez Magro (1979) (<», McCarthy Howell & Low (1980) (A), 
Mariotti et al. (1983) (o), Dyck et al. (1984) (+), Dyck et al. 
(1987) (*), Benson, Turner & Dyck (1989) (®) and Danchi et al. 
(1990), (1994) (□). Phases and angular sizes of the dust conden- 
sation zone are listed in Table 2. 



Table 2. Angular size of the dust condensation zone at maximum 
and minimum light determined from the available visibility obser- 
vations. References for the observations are listed in the caption 
of Fig. 2. 
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<J)(a) 






(^m) 




(arcscc) 


(per cent) 


2.2 


max 


0.40 


29 


2.2 


min 


0.29 


17 


5 


min 


0.32 


8 


10 


max 


0.45 


20 


10 


min 


0.30 


14 



(a) Phase of the light curve at which visibility observations were 
made at the wavelength listed in the first column 

(b) Angular size of the dust condensation zone determined from 
model fit to the visibility observations 

(c) The percentage difference between the angles listed in column 
(b) and those determined from equation (3) and the best fit to 
the SED (0.56 arcsec for max and 0.35 arcsec for min). 

from the spectral shape we fit our models to visibility ob- 
servations at minimum and maximum light obtained at 2.2, 
5 and 10 vim. In these fits, 9i is taken as a free parameter, 
providing an independent estimate for it. Fig. 2 shows com- 
parison of model results with observations. Since the data 
at 2.2 \x.m are spatially asymmetric, the plotted results are 
spatially averaged. Phases and values of 6\ are summarized 
in Table 2. The independent fits for di from the visibility 
and the SED agree within 15-20 per cent on average. It can 
be estimated that, within 20 per cent, the angular size of the 
dust condensation point varies between 0.3 and 0.5 arcsec. 

Previous models have never achieved simultaneous 
agreement for both the SED and spatially resolved obser- 
vations at short wavelengths (e.g. Martin & Rogers 1987). 
In all these models, optical depths that fitted the SED pro- 
duced a 2.2-ii-m visibility too large at g ^ 1 arcsec"^. This 
problem is directly related to the flux shortage of the mod- 
els at these wavelengths and is another manifestation of 
the need for large grains. We have shown in IE96 that the 
value of the visibility when it levels off at large q is simply 
exp(— Tsca) for A ^ 3 ^.m, where Tsca is the scattering optical 
depth. Therefore, Tsca must be increased to reduce the vis- 
ibility. Also, because Tabs is fixed from the spectral shape, 
this increase translates to an increased albedo, implying the 
presence of grains with sizes of ~ 0.2 (i.m. This independent 
estimate of the grain sizes provides further support for the 
one obtained from the spectral shape. 

Recently, Danchi et al. (1994) obtained visibility curves 
for IRC-f 10216 at 11 \ira close to maximum and minimum 



Figure 3. Single-scan (E-W) imaging of IRC-l-10216 at 10 [am. 
The thick solid line in the top panel is the observations of 
Bloemhof et al. (1988). Superimposed on it is our model result 
drawn as a dashed line, hardly distinguishable from the obser- 
vations. It is obtained by a two-dimensional convolution of the 
surface brightness for d\ = 0.35 arcsec (innermost thin solid line) 
with the point-spread function (PSF, dot-dashed line; all profiles 
are normalized to unity at the peak). In the bottom panel, the 
dashed line is the surface brightness deduced by Bloemhof et al. 
by one-dimensional deconvolution of the observed profile with the 
PSF. The dot-dashed line is the one-dimensional convolution of 
our model result with the PSF. The thick solid line is the two- 
dimensional convolution of our model result for 6i = 0.35 arcsec 
with the PSF (the same as the dashed line in the top panel). 
The two thin solid lines below and above this curve correspond 
to curves derived analogously from models with d\ = 0.30 and 
0.40 arcsec, respectively. 



light. From these observations they find 6\ ~ 0.1-0.2 arcsec. 
Til = 1.24 and T\ = 1360 K. These results, determined by 
fitting visibility curves at the single wavelength 11 [xm, differ 
from ours; by comparison, from the spectral shape we find 
Til = 0.3 — 0.4 and T\ — 750 K. It is important to note that 
the Danchi et al. data, displayed as solid and open squares in 
the bottom panel of Fig. 2, are properly fitted by our models 
together with all other data. By contrast, the Danchi et al. 
fits rely on a limited data set, confined to a single wavelength 
and visibilities < 0.4, which does not sufficiently constrain 
the model parameters (IE96). Indeed, the values of th and 
Ti deduced by Danchi et al. cannot produce a simultaneous 
fit to the visibilities and SED at wavelengths shorter than 
5-6 [im. 

Bloemhof et al. (1988) obtained a single-scan image of 
IRC-l-10216 at 10 ^im close to minimum light (phase ~ 0.4). 
We computed the profile expected in those observations from 
the model surface brightness determined for this phase from 
the spectral shape. The top panel in Fig. 3 shows the com- 
parison between the observed profile (outermost thick solid 
line) and our model result (dashed line, overlapping the ob- 
servations; all displayed profiles are normalized to unity at 
their peaks) The innermost thin solid line is the surface 
brightness distribution obtained from our model with Oi — 
0.35 arcsec. The central peak corresponds to the stellar con- 
tribution and the features at relative RA ±0i /2 to the dust 
formation zone. The model result is obtained from a two- 
dimensional convolution of this profile with the point-spread 
function (PSF), shown with the dot-dashed line. This fit to 
the observed image provides independent determination of 
9i around minimum light, in agreement with the previous 
two. 

Our model surface brightness is considerably different 
from the one deduced by Bloemhof et al., displayed with 
the dashed thick line in the bottom panel. In addition to 
the contribution of a central component with a width of ~ 
0.4 arcsec, this profile also requires a more extended com- 
ponent with a width of ~ 2 arcsec. Bloemhof et al. obtained 
this distribution from a one- rather than a two-dimensional 
deconvolution of their observed profile with the PSF. We 



t To remove a slight asymmetry in the observations, all profiles 
are symmetrized east— west. 
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have verified that the results of the commonly used one- 
dimensional convolution with the PSF are indistinguishable 
from those of the proper two-dimensional convolution for 
centrally peaked surface brightness. However, in the case of 
extended structures whose peaks do not coincide with the 
centre of symmetry (e.g., a ring), the two procedures pro- 
duce different results. A one-dimensional convolution of our 
model surface brightness with the PSF produces the thin 
dot-dashed line in the bottom panel, considerably different 
from the result of the proper two-dimensional convolution, 
repeated in the bottom panel as the thick solid line. 

The possible existence of an extended ~ 2 arcsec com- 
ponent was first conjectured from lunar occultation observa- 
tions by Toombs et al. (1972) because they could not prop- 
erly fit their results with a single, sharp-edged disc with a 
diameter of ~ 0.4 arcsec. Thus they invoked an additional, 
larger sharp-edged disc. Recently, Sloan & Egan (1995) also 
modelled their 10-fi.m observations, obtained with 0.9 x 2.0 
arcsec^ resolution, in terms of single-temperature compo- 
nents, a procedure that produced two dust shells. One shell 
ranges in diameter from 0.055 to 0.67 arcsec with a single 
temperature of 340 K, the other from 1.5 all the way to 
5.2 arcsec with a single temperature of 240 K. However, our 
model, which has no sharp edges or discrete temperature 
components, properly explains both the Toombs et al. and 
Sloan & Egan observations because it is extended. There- 
fore, our modelling does not support the existence of the 
conjectured 2 arcsec component or any other discrete shell. 
We find no need to augment the steady-state outfiow with 
any additional components. 

Direct imaging is a most sensitive method for determin- 
ing the dust condensation radius of optically thin envelopes. 
In the bottom panel of Fig. 3 we display in thin solid lines 
the imaging results expected with the PSF of Bloemhof et al. 
when 9i is varied by only ± 0.05 arcsec. Since the expected 
variation of 6\ between minimum and maximum light is con- 
siderably larger, measurements of this variation during the 
stellar cycle can provide an important check of our models. 



4 DYNAMICAL PROPERTIES 

The quantities determined from IR observations can be used 
to constrain dynamical properties of the IRC+10216 out- 
fiow. Momentum flux conservation relates the mass-loss rate 
M = M-5 X 10~^ Mq yr~^ and Uo, the terminal outflow ve- 
locity in km s^^, to radiative properties via 



M_5«e = 20rFL4(l - 



(4) 



This relation is valid in steady-state for the time-averages of 
tf and the stellar luminosity L» = 1/4 x 10* Lq , where F (oc 
Mt/Li) is the gravitational correction (IE95). Our model 
calculations give F"^ = 0.13 for IRC-l-10216 if the stellar 
mass is M* = 1 M0 and L4 = 1.5; uncertainties in M* and 
L4 are of minor importance since F~^ is so small. From our 
models we flnd that rp varies from 2 at maximum light to 2.4 
at minimum (see Table 1). With a terminal outflow velocity 
of 15 km s~^ (Morris, Lucas & Omont 1985; Zuckerman & 
Dyck 1986), this gives 



where L4 refers to the luminosity at maximum light, assum- 
ing a bolometric amplitude of 1 mag. 

The dust mass-loss rate M is directly related to the 
optical depth via 



3Af 



(6) 



Q(a) 



n,j(a) da. 



Here pa is the dust solid density (1.85 g cm~^. Rouleau 
& Martin 1991); Q{a) is the extinction efficiency of grains 
with radius a; rid is the dust density profile normalized to 
unity at ri; and V\ is the outfiow velocity at ri, assumed 
to be 1 km s~^. This velocity corresponds to the velocity 
at the sonic point (Deguchi 1980) and introduces the prin- 
cipal uncertainty in determining M"^ . From this result and 
expressions for ri listed in IE96, our model calculations of 
IRC-H0216 produce 



M_ 



■ 4.2a/I^ 



(7) 



where M.g = M /lO"* M© yr"^; the variation with y/Tl 
refiects the dependence of ri on luminosity (cf. Section 3). 
The ratio of the mass-loss rates of the entire envelope and 
the dust component produces the average gas-to-dust ratio 



1000 



M_ 



= 500Vl^ 



(8) 



In all these relations, the luminosity remains unknown. 
Since the bolometric fiux is 2.1 x 10~* W m~^, the luminosity 
obeys 



L4 = 1.5 Dfs 



(9) 



where D — Diso x 150 pc is the distance to the star. In 
terms of this parametrization, the various quantities listed 
above are 



M_5 = 3.1D?5o, M_s^5.2Di 



= 610i3i5o .(10) 



These results are in good agreement with independent esti- 
mates if we take D — 150 pc. Estimates for M_5 based on 
CO data range from 1-2 (Jura 1994, and references therein) 

to 4.7 (Kwan & Webster 1993). Values derived for M^g vary 
from 5.6 (Grifiin 1990) to 9.8 (Kastner 1992). Gas-to-dust 
ratios for carbonaceous winds in late-type stars range from 
260 (Volk, Kwok & LangiU 1992) to 670 (Knapp 1985). 

The most-often quoted distance to IRC-f-10216 is 290 
pc, derived by Herbig & Zappala (1970) for an assumed 
luminosity of L4 — 5.5. This distance leads to unrealisti- 
cally high mass-loss rates and gas-to-dust ratios. A closer 
distance of 100-150 pc has been proposed by Zuckerman, 
Dyck & Claussen (1986) and Kastner (1992), and our results 
support these suggestions. Based on theoretical considera- 
tions, Martin & Rogers (1987) pointed out that L4 ^ 1.2. 
Thus, IRC-f 10216 is probably not closer than 130 pc, and 
the value of 150 pc adopted by Jura (1994) and in this work 
is likely to be close to the true distance. With this distance 
of 150 pc, the radius of the dust condensation zone varies 
from ~ 3 X 10^* cm at minimum light to ~ 5 x 10^* cm at 
maximum. 



M-5 = 2.IL4 



(5) 
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5 DISCUSSION 

Our model for IRC+10216 provides a description of the IR 
observations based on a self-consistent treatment of the dy- 
namics and radiative transfer. The model provides simulta- 
neous agreement for both the SED and all high-resolution 
observations of this star. Particularly encouraging is the 
agreement between the three independent determinations of 
the angular size of the dust condensation point (~ 0.3 arcsec 
at minimum light). Our results demonstrate the advantage 
of following the two-step modelling procedure outlined in the 
Introduction. The most appropriate approach is first to con- 
strain the overall optical depth by the spectral shape, then 
use visibility data to determine the size of the envelope. 

Because of the close agreement obtained for such a va- 
riety of independent observations, we do not expect major 
changes in the parameters determined here. The model could 
still be improved by considering the asymmetry of the en- 
velope. As mentioned in Section 2, although symmetric at 
wavelengths longer than ~ 3-4 jxm, the observed surface 
brightness is slightly asymmetric at short wavelengths. Such 
a dual appearance can be understood in terms of the ba- 
sic physical processes that control the IR radiation at the 
different spectral regions. At short wavelengths, scattering 
dominates the observed radiation. Since scattering can be 
expected to map the entire envelope, scattered radiation 
should reflect the density distribution, displaying any asym- 
metry in it. As long as the elongation is not severe, it is 
reflected only in the shape of the image, not in the flux. 
On the other hand, radiation at longer wavelengths is dom- 
inated by dust emission, predominantly controlled by the 
dust temperature distribution. For slightly elongated den- 
sity distributions, the dust temperature distribution can still 
be spherically symmetric to a good degree of approximation 
because it is mostly controlled by the distance from the cen- 
tral star. This explains the close agreement of our spherically 
symmetric model with the data. A slightly enhanced mass- 
loss rate in the equatorial plane, as proposed by Ridgway 

6 Keady (1988) and supported by Kastner & Weintraub 
(1994), can be accommodated without a significant effect on 
our results. On the other hand, our model cannot describe 
the emission from a bipolar nebula, the geometry suggested 
for IRC-l-10216 by Dyck et al. (1987), and its success indi- 
cates that such a drastic departure from spherical symmetry 
may not be necessary to explain the observations. 

We resolve the difhculties encountered in previous stud- 
ies at short wavelengths by including large grains. The short- 
wavelength behaviour of both the SED and the visibility 
shows that grains as large as ~ 0.15-0.2 \im are present. 
Jura (1994) recently discussed the grain size distribution 
for IRC-l-10216, based on polarization in the K band and 
shielding of circumstellar molecules against destruction by 
interstellar UV radiation. He finds that grains as large as 
~ 0.1 \im exist in the outer envelope (more than 15 arcsec 
from the star), in good agreement with the sizes obtained 
here. Since our analysis of grain sizes is primarily affected 
by the inner regions of the envelope while Jura's results ap- 
ply to the outer regions, we conclude that the grain sizes 
do not change signiflcantly through the envelope. The ef- 
fects of sputtering, grain growth, etc., do not seem to be too 
important after the initial dust formation. 
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